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The semi-annihilation reaction takes the schematic form V'/V^/ ^ Wk0^ where i//, are stable dark 
matter particles and is an unstable state. Such reactions are allowed when dark matter is sta- 
bilized by a larger symmetry than just Z2. They lead to non-trivial dark matter dynamics in the 
early universe, and the thermal production of the relic particles can be completely controlled by 
semi-annihilations. This process might also take place today in the Milky Way, enriching the 
(semi-)annihilation final state spectrum observed in indirect detection experiments. 
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1. Introducing semi-annihilation 

The evidence for dark matter is one of the best motivations for physics beyond the Standard 
Model (SM) [1]. A well-motivated dark matter candidate are so-called Weakly Interacting Massive 
Particles (WIMPs), whose production in the early universe is thermal [2, 3]. The thermal freeze-out 
is very predictive, thus it is important to know how to correctly compute the relic density, especially 
when the relic computation cannot be reduced to the standard case (see e.g. [4]). 

The thermal abundance can be dramatically affected by the presence of the "semi-annihilation" 
reaction. If the dark matter is composed of more than one stable component y/,, it takes the form 

YiYj^Yk^, (1.1) 

where is a SM state or a new particle which decays to the SM. We see that unlike ordinary 
annihilation where the total dark matter number changes by two units, in semi-annihilation the 
total dark matter number changes by only one unit, as shown in figure 1 . As long as the triangle 
inequality nik < mj + mj is satisfied (as well as its crossed versions), the semi-annihilation reaction 
— Yk'P can take place without making any relic particle unstable. Such reactions take place 
once we relax some assumptions about the symmetry structure of WIMP interactions. For example 
the dark matter can be stabilized by a Z3 symmetry [5] , composed of non-Abelian gauge bosons 
[6], or more generally stabilized by "baryon" and/or "flavor" symmetries, as in QCD-like theories. 
We study two expUcit toy examples of such models [7]. We find that to correctly compute the relic 
abundance, semi-annihilation must be included, and it can dominate over standard annihilation. We 
also explore the effects of semi-annihilations on indirect detection experiments. 

We introduce semi-annihilation with a simple case. We assume dark matter to be composed 
of a complex scalar x, which is stabilized by a Z3 symmetry. The x particle interacts with a real 
scalar 0, which eventually decays to SM states. The x particles are kept in thermal equilibrium in 
the early universe by the reactions XX ~^ 00 ^^'^ XX ~^ X^- While we focus on the case of a field 
(p "portal" to the SM [8], in more general dark matter scenarios (j) could be a SM field itself [9]. 

Single-component dark matter with a Z3 symmetry is not representative of models in which 
semi-annihilation is relevant. The second model we consider is a minimal multi-component dark 
matter model having semi-annihilations among its allowed reactions.' In this toy system, dark 
matter is composed of two stable components, a complex scalar field X and a vector-like fermion 
b and f^. We impose a.U{l) global symmetry under which b, b'^, and x have charges +1, —1, and 
—2, respectively. We also introduce a real scalar- portal ^ field, which we assume to be in thermal 
equilibrium in our calculation and decays to SM states. Fermion number conservation always 
guarantees that b is stable, whereas stability for x requires the triangle inequality rule rii), < Inih 
(isosceles triangle). The model and its symmetries are summarized in table 1. 

The remainder is structured as follows. In Section 2 we present a study of the relic density 
in the two toy models we consider. In Section 3 we explore the effects of semi-annihilations on 
indirect detection experiments, and conclude in Section 4. 

'a more realistic construction, where a supersymmetric gauge thieory with N f = + I is considered, can be found 
in [7]. Here, we prefer to focus on a simpler model avoiding unnecessary complications caused by having superpartners. 
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Figure 1: Reactions that keep the dark matter particles in thermal equilibrium in the early universe: (a) 
annihilation; (b) semi-annihilation. The field decays to SM states. 
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Table 1: Field content and symmetries of a minimal multi-component model with semi-annihilation. 

2. Implications for the relic density 

The number density evolution of % in the model with Z3 symmetry is described by ^ 
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(2.1) 



In the {ov)j^j^^j^^ = limit we recover the standard case. The above Boltzmann equation can be 
solved semi-analytically in analogy with the method in [3], by solving the equation in two regimes, 
early and late times, and then matching the two solutions at the freeze-out point. The freeze-out 
temperature Xf, where we define x = m^/T, is found by solving the equation 



Xf = log 
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(2.2) 



where c = Vl — 1 for s-wave amplitudes [3]. The first term on the right-hand side of (2.2) is 
the solution we usually have in the Lee-Weinberg scenario. Thus we see that the effect of semi- 
annihilation is to shift the freeze-out temperature by only a small logarithmic amount. The mass 
density of the x particle today is given by 

1.07 X lO'^GeV-i 



nyh-" = 2 X 



J{Xf) 



dx 
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1 , , 
+ 2\'^'^ixx^x<^ 



(2.3) 



In single component dark matter models, the conventional assumption is that we have only an- 
nihilation (namely = 0), thus once we fix the dark matter mass m-^^ the cross section 

^We correct the Boltzmann equation in [7] with the 1/2 factor in the semi-annihilation term. We thank Brian Batell. 
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(a) (b) (c) 

Figure 2: Y vs. x evolution in the bbx model for the b (blue lines) and the % (green lines). We define 
Yi = ni/s, where i = b,X s is the entropy density. We plot the equilibrium values (dashed lines) and the 
numerical solutions (solid lines). We fix m^^ — 0.8 TeV, nij, — ITeV, a —2 and j3 = 0.01. The values of fc, 
e and i^Divi/i^ are shown in each plot. Here, x particles would be over produced (a), but because of phase 
space suppression K (b ^ x conversion) is ineffective (b). However, e (semi-annihilation) is never phase 
space suppressed and can restore the desired relic density (c). 

{av)),2^(j)il, is uniquely determined by the relic density measurement. However, if we allow also 
semi-annihilations, a wider region of parameter space gives the correct relic density, thus the ther- 
mal production of the relic particles can also be completely controlled by semi-annihilations. 

In models which involve more than one stable dark matter component a semi-analytical so- 
lution for the relic density is simply not possible, and one must resort to numerically solving the 
complete set of coupled Boltzmann equations. In particular, the simplifying assumptions made 
when analyzing co-annihilating models [4] are not applicable for semi-annihilation [7] . As a con- 
sequence, the dark matter dynamics in the presence of semi-annihilation are far more varied than 
in decoupled multi-component scenarios. 

In our second example model, the relic abundance of b and x is described by a system of 
two coupled Boltzmann equations. We can now write down and solve numerically the Boltzmann 
equations for the relic abundance of b and x and see how semi-annihilations affect the dark matter 
relic density. Once we fix the dark matter masses and limit our consideration to 5-wave annihilation, 
there are four free parameters, which we choose to be the following s-wa.ve matrix elements 

•^bb^^^ = •^XX^'^i' ~ ^' -^bb^xX ~ ^^xb^^b ~ ^' 

The first two amplitudes come from the standard annihilation to light particles in the final state, the 
third from the process where two dark matter particles of one species are entirely converted to two 
dark matter particles of the other species. This process is phase space suppressed for typical kinetic 
energies at freeze-out. The last corresponds to the semi-annihilation process, and we are interested 
to see how the relic abundance changes with the parameter e. 

In our numerical study we fix the dark matter masses to be m;^ = 0.8 TeV and mi, = ITeV 
for concreteness, and we vary the matrix elements amplitudes. We fix the values of a and j8 
(annihilation amplitudes) and study how the relic density is affected once we turn on the parameters 
KQy^ X conversion) or e (semi-annihilation). The main difference between semi-annihilation and 
species conversion is that the former is never phase space suppressed. This is beautifully illustrated 
in the case when pure annihilations would give a lot of x and very few b, as shown in figure 2(a) 
for a = 2 and j3 = 0.01. The semi-annihilation process x^ — )• b<^ is still able to destroy x particles 
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Figure 3: Two different contributions to indirect detection spectra: (a) annihilation; (b) semi-annihilation. 
The field (psM is a portal field, but in more general dark matter scenarios could be a SM field itself. 

and bring the relic density back to the observed value, figure 2(c) for £ = 5. However, the reaction 
XX ~^ is powerless, even if we badly break perturbation theory for = 4;r as in figure 2(b). 
This reaction is forbidden at zero kinetic energy, and at freeze-out the thermal energy is very small 
compared with the mass splitting between the two particles. Thus, semi-annihilation is a truly 
unique species changing interaction that affects early universe cosmology. 

3. Implications for indirect detection 

We now explore the implications of semi-annihilations for dark matter detection experiments. 
Since semi-annihilation does not give any new contributions to direct detection rates, we focus on 
indirect detection only. There has been much recent interest on indirect detection of dark matter, 
motivated by observations suggesting the presence of a new primary source of galactic electrons 
and positrons [10]. These anomalies might be evidence for dark matter annihilation or decay in our 
galaxy. 

In general, semi-annihilation does not lead to the "boost factor" necessary to explain the mag- 
nitude of these candidate dark matter indirect signals. However the presence of such a reaction 
can still have considerable implications for indirect searches. Semi-annihilations are additional 
channels to produce light particles in dark matter interactions in our galaxy, as shown in figure 
3, which have no analog in standard multi-component scenarios. Thus the predicted spectrum is 
enriched with respect to the standard scenarios. In addition, our examples models can account for 
the fact that dark matter (semi-)annihilates preferably into leptons; if we assume a portal (j) field 
mass < 2GeV the decay of (p to antiprotons is kinematically forbidden [8]. 

We present here the results for the bbx model. ^ In this case, there are four reactions to pro- 
duce (/> particles in the final state, the ordinary annihilations of the two components and two semi- 
annihilations. Thus the spectrum is this multi-component case is quite rich, since the four reactions 
yeld four monochromatic <p lines. We present numerical results for the (/> spectrum for the same 
values of the dark matter masses chosen in Section 2, namely = O.STeV and nih = ITeV. To 
simplify the following discussion we perform our analysis for K" = 0. Two example spectra are 

^For the spectra in the Z3 model see reference [7]. 
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Figure 4: Spectra of (p particles in the bbx model for thermal production mostly controlled by: (a) annihi- 
lations, (b) semi-annihilations. The values of the matrix element amplitudes are shown in each plot. 



shown in figure 4, where two opposite cases are considered. In the first case, figure 4(a), the 
semi-annihilation contribution is taken to be very small, thus the pure annihilation lines dominate. 
When the thermal production is mostly controlled by semi-annihilations, as in figure 4(b), the 
semi-annihilation lines overwhelm the standard contributions. 

One expects the indirect detection spectra in models with more than two stable components 
to have an even richer structure of lines. Of course, once convolved with the decays, the 
spectrum is less distinct [7]. More detailed study is necessary to know whether such a spectrum 
can be measured in 7 rays telescopes. 



4. Conclusions 

We relaxed some assumptions about the symmetry structure of WIMP interactions by allowing 
relic particles to "semi-annihilate". We studied two explicit examples where semi-annihilation is 
present: a single species dark matter model with a Z3 symmetry and a multiple species dark matter 
model with "baryon" and "flavor" symmetries. We saw that the semi-annihilation process can 
have a considerable effect on the dark matter relic abundance, and the dark matter dynamics in the 
presence of semi-annihilation are far more varied than in decoupled multi-component scenarios. 

The inclusion of semi-annihilation has interesting implications for indirect detection experi- 
ments. The overall integrated flux is not very affected by semi-annihilation. However, the final 
state spectrum in semi-annihilating models is far richer than in standards scenarios because of the 
differing kinematics between semi-annihilation and ordinary annihilation. 

In the language of [4], semi-annihilation is in some ways the "fourth exception" in the calcu- 
lation of dai^k matter thermal relic abundances. We find it intriguing that unlike the traditional three 
exceptions (co-annihilation, annihilation below a mass threshold, and annihilation near a pole), this 
fourth exception not only affects dark matter interactions in the early universe, but also leaves an 
imprint today via the indirect detection spectrum. We expect that there are a wide variety of multi- 
component dark matter models with species changing interactions, motivating further studies of 
the semi-annihilation process. 
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